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ground of a non-Abelian Weizsacker- Williams field. These corrections are valid to 
all orders in the effective coupling a s fi, where [i 2 denotes the average valence quark 
color charge squared per unit area. We find ln(l/x) logarithmic corrections to 
the classical gluon distribution function. The one loop corrections to the classical 
Weizsacker- Williams field do not contribute to these singular terms in the distribu- 
tion function. Their effect is to cause the running of a s . 



1 Introduction 



In a recent work [jj]] we computed the Green's function in light cone gauge A + = for 
the small fluctuations about a background Weizsacker-Williams gluon field. This 
background field is generated by the valence quarks in a large A nucleus. For small 
x partons, these valence quarks constitute a static and well localized source of color 
fields ||]. The average color charge squared per unit area of the valence quarks 
is denoted by and it is of the order of A 1 / 3 fm -2 . The quantity fi is the only 
dimensionful parameter in the theory and as a result, the coupling constant will run 
as a function of it. 

Previously, two of us (McL.-V.) have argued that we can compute the gluon 
distribution function from the light cone gauge Green's function ||. In the present 
paper we use that Green's function to compute the corrections induced by quantum 
fluctuations on the Weizsacker-Williams distribution function. 

There are several reasons for this computation to be of interest. From the 
practical point of view, we hope that an analysis which provides us with a better 
understanding of the initial conditions for the evolution of partons in the collision of 
heavy nuclei to form a quark-gluon plasma will help establish a firm foundation Q 
for partonic cascade models simulating such collisions ||. For an alternative ap- 
proach to the problem of initial conditions in heavy ion collisions, see reference ||. 

From the theoretical point of view, we hope to understand the small x behavior 
of a nucleus starting from a QCD based approach. Let us recall that it is believed 
that in the small x region, the gluon distribution function computed perturbatively 
including leading logarithmic contributions in ln(l/a;) for a single nucleon behaves 
like § 

dN 1 

dx x 1+Cas 

The steep rise of the gluon distribution function for small x is sometimes referred to 
as the Lipatov enhancement. It is obtained by solving the Balitsky-Fadin-Kuraev- 
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Lipatov (BFKL) kernel for the t -channel exchange of a perturbative pomeron 
This behavior is also exhibited in a hadron where the large x part of the hadron 
wavefunction is taken to be a heavy quark-antiquark state. By applying Hamilto- 
nian perturbation theory to this state, it is possible to reproduce the kernel of the 
BFKL equation H for the emission of a large number of soft gluons. 

However, the rapid rise of the gluon distribution function with smaller values 
of x is in conflict with unitarity when considering the hadron scattering total cross 
section at asymptotically high energies 0. Physically, this violation of unitarity 
can be understood to result from ignoring effects that arise due to the large density 
of partons at very small values of x |10(| . When the density of partons is so large 
that neighboring partons overlap, the t-channel picture of an independent parton 
cascade in x breaks down. The former signals that at very small values of x, the 
picture in which the partons do not interact with each other has to be modified in 



order to comply with the Froissart bound [11 1 on the growth of cross sections at 
asymptotically high energies. Although some work has been done in recent years 
to include these "higher twist" effects in describing parton evolution at high densi- 
ties |ll], 12, 13], more still remains to be done in devising a quantitative mechanism 
to limit this growth. 

The regime of high parton density is the screening regime. This screening is 
presumably responsible for the shadowing phenomena observed in deep inelastic 
scattering experiments off nuclei at small x. It can be addressed as a collective or 
many body effect. It is precisely this many body problem of parton interactions that 
we seek to address in our work. As outlined in references [|| fli!] ], our formalism 
provides us, by means of a novel weak coupling approach, with a technique to solve 
the many body problem of wee parton distributions in a large nucleus. 

In this work, we will focus on one of the theoretical aspects of the problem -the 
nature of the small x terms in the gluon distribution function of the Weizsacker- 
Williams fields generated by the valence quarks in the infinite momentum frame. We 
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do this by computing a formula for the gluon distribution function which includes 
the effects of all orders in the parameter a s jJL. Working in the weak coupling regime 
a s fi <C kt, we extract from this formula an expression for the distribution function 
in perturbation theory to second order in a s . Although we approach the problem 
in weak coupling, we will show that after the corrections are considered, the series 
expansion parameter becomes a s ln(l/x) which may be large for small x's. We 
will argue that this forces us to devise a method to isolate and sum up the leading 
contributions to all orders in that effective expansion parameter in order to compute 
the modification to the zeroth order 1/x distribution function. 

The outline of this work is as follows: In section 2, we briefly review the basic 
aspects of the model in which we treat the nuclear valence quarks as static sources 
of color charge as seen by small x partons. We also go briefly through the formalism 
that allows us to compute the Green's function for the small fluctuations equation 
and finish the section by writing the formula for this Green's function. In section 
3, we use this Green's function and exploit its relation to the gluon density in order 
to compute the gluon distribution of small x gluons. We also derive from this result 
a formula for the leading small x terms of the distribution function in perturbation 
theory valid to second order in a s . In section 4, we compute the corrections to 
the background field induced by the fluctuation field and show that the only effect 
that this introduces can be absorbed into the renormalization of the background 
field (which is related to the renormalization of the coupling constant to one loop). 
Details of the calculations in section 3 and section 4 are discussed in appendix I and 
appendix II respectively. Finally we summarize our results in section 5. 

2 The model 

In QCD, a hadron is a cloud of virtual particles with a rather complicated struc- 
ture. The picture gets simpler when we consider the hadron as a quantum system 
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composed of quasi real particles (partons) with lifetimes much larger than the char- 
acteristic interaction times. This can be done in a reference frame where the hadron 



has a large momentum [10]. Partons with large lifetimes can produce new partons 
carrying smaller fractions x of the initial hadron's momentum. The small x partons 
will therefore densely populate the hadron and see the rest of it with its longitudinal 
dimension Lorentz contracted to a thin disk. In our model, we look at the small 
x partons in a large A nucleus (x <C A -1 ' 3 ) where the high parton density allows 
us to use weak coupling techniques. The rest of the nucleus consists of the valence 
quarks which carry most of the nuclear momentum. They are described as a static 
(recoilless) source of color charges, in a reference frame in which they move with 
the speed of light (infinite momentum frame) [0]. 



The problem is well suited to be described using light cone variables [15| 



y°,y 3 -^y ± = (y°±y 3 )/V2 (1) 

In order to compute ground state properties of the wee partons, we define a 
partition function for the system. This partition function includes the sum over 
a large number of color configurations. To simplify the problem, we resolve the 
transverse space direction into cells which contain a large number of valence quarks, 
or equivalently, a large number of color charges. This allows us to treat the sum 
over color configurations classically Q. To write the average over the color charges, 
we introduce a Gaussian weight by inserting into the path integral representation 
of the partition function the term 

exp(--^ fd 2 x t p 2 (x)\ , (2) 



v 2p 2 , 

where p is the color charge density (per unit area) and the parameter p? is the 
average color charge density squared (per unit area) in units of the coupling constant 
g. The introduction of the partition function, where we average over the sources of 
external charge, allows us to formulate the theory as a many body problem with 
modified propagators and vertices. 



We treat the system perturbatively and the first step is to solve the classical 
equations of motion 

D^ = g r, r = 5 +u p(x + ,x t )5(x-) 1 (3) 

for which (working in the light cone gauge A_ = —A + = 0) there exists a solution 
with A + = —A~ = 0. We require 

Ai(x) = 9(x-)ai(xt) , (4) 

(hereafter, latin indices refer to transverse variables) and furthermore F lJ = 0. The 
latter condition implies that a(xt) is a pure gauge transform of the vacuum 14], 



T-a i = -U(x t )V i rf(x t ), (5) 
9 

where U(xt) is a 577(3) local gauge transformation whose spatial dependence is 
only on the two-dimensional transverse space. It is subject to the physical gauge 
condition 

V-(u(x t )VU\x t ))=-ig 2 p. (6) 

The x + (light cone time) dependence of the charge density is a consequence of the 
extended current conservation law 

D„ J M = 0. (7) 

The integration over the sources p in equation (|2|) may be written as 

(8) 



[dU] exp f-^jl* • U ] ViU ^ 



where we have ignored the Faddeev-Popov determinant. Note that the effective 
coupling constant for this theory is g 2 p so that the expansion parameter becomes 

OLsP/Pt- 
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The Green's function can be computed from the relation 

GfM =1^-1 ^^-p 2 )WUx){AP){{x') (9) 



X-ieJ (2vr) 4 

and the gluon distribution function can be computed from the Green's function by 
the relation 

fJAT <2J?+ 

ft=^EC(^,fc;^fc). (io) 

To relate the Green's functions to the distribution function by the above relation, 
the former must be averaged over the external sources of color charge. 

Indeed, the distribution function, to all orders, is related by the above expression 
to the fully connected two point Green's function. This Green's function is given 
by the relation 

((AA)) p = ((A cl )(A cl ) + (A q A q )) p . (11) 

In the above, (Ad) is the expectation value of the classical field to all orders in U. 
It can be expanded as 

(A cl )=A^+A^ + , (12) 

where A^' is the solution discussed in Eqs. (4)-(5). The one loop correction to the 
classical field, A^\ is computed in section 4 of this paper. The term (A q A q ) above 
is the small fluctuations Green's function computed to each order in the classical 
field. The symbol < • • • > p indicates that we have to average over the external 
sources of color charge with the Gaussian weight described above. 

From the above, it is clear that the zeroth order contribution is (A^ A^ 1 ) p . 
This contribution is the QCD analog of the well known Weizsacker-Williams dis- 
tribution in classical electrodynamics. The general form of the solution is given in 



reference [14]. In the range of momenta ct s /i <C kt <C the zeroth order solution 
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yields a distribution function that, written in terms of x = k + / P + , with P + the 
nuclear longitudinal light cone momentum, looks like 

1 dN _ a sf , 2 (N 2 -l) 1 
ttR 2 dxcPkt 7r 2 xk 2 

The above is the well known Weizsacker- Williams distribution scaled by fx 2 ~ A 1 / 3 
fm~ 2 . 

With this formalism at hand, we can proceed to compute the next order con- 
tribution to the gluon distribution function. Our strategy is to compute the small 
fluctuations correction to the classical equation of motion. Writing the field in terms 
of its background and fluctuation parts 

A tl {x) = B^{x) + b^(x), (14) 

we are able to express the equation obeyed by 6 M as ]l]] 

V " D^(B)D u (B)^j b u - 2F^b v = 0. (15) 

where B is the background field which according to equation (||) is non-vanishing 
only for its transverse components. D fl {B) is the covariant derivative with B as 
the gauge field (notice that = d^). As discussed in reference Q, the set of 
equations (^) can be unambiguously solved in the gauge A~ = 0, and by means of 
equation @ we can compute the Green's function for the fluctuation fields in this 
gauge. To obtain the Green's function in the gauge A + = (light cone gauge), we 
perform a gauge transformation on the Green's function in the A~ = gauge and 
obtain finally (in the matrix representation) 



aim r d A n e^-y) 

Gf^(x,y)= f dpe 



| (2c ip+ ( x ~ y } — e~ ip+y —e ip+x 

(2-7r) 4 p 2 — ie 1 L 13 p~p + 

e(-x-)e(-y-)rfT^' + e(x-)9(y-)F^(x t )Ft' P '(yt) 

+ 0(-x-)8(y-) [ j-%d 2 z t e *(« + -P + )wV(pt-«)(i«-*0 
J (2tt) 



6 + P^L(jp + *- 
p p+ 



+ ML, e - iq + y - PiQjPyqt ( e i P +x- _ 1)(e -^ + y- _ V 

P -q+ (p~p + )(p~q + ) 

+ 9{x-)9(-y-) f-f^d 2 z t e i(pt~<it)(z t -x t ) e -i( q +- P +)x- 

J (27T)^ 



QiPjPt ■ qt 



Si . + mL( e -^ + y-- 1) + m_^- _ 1} 



p p 



(p p+){p q+) 



e -vrv -l)( e iq x - I) 



(16) 



r,2 „2 



where q + = p + + 9i 2p Z' i and 



FW(x t ) = {u{x t )r a U\x t )Y . (17) 



The above result was derived in reference |M. The regularization of the poles in 
p~ and p + is a rather subtle issue and was discussed at some length in the above 
mentioned paper. 

The Green's function Gij(x,y) contains the physical degrees of freedom that 
we need to relate to the gluon density. In the following section, the above small 
fluctuations propagator will be used to compute corrections to the Weizsacker- 
Williams distribution function. We will show that logarithmic corrections, both in 
x and kt, arise from here. In section 4, we discuss the corrections that the small 
fluctuations induce on the classical field to one loop. These corrections, in contrast, 
provide no such logarithmic terms to the distribution function and their only effect 
is to renormalize the coupling constant and the background field. 

In appendix I, we compute the Fourier transform of the Green's function. This 
result is useful for the computation of the distribution functions performed in the 
next two sections. 
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3 The gluon distribution function 

In this section, we will show how one can obtain the gluon distribution function 
from the Green's function (equation [l(]), corresponding to the small fluctuations 
about the Weizsacker-Williams background color field. We will be concerned with 
the structure of the leading terms for small x and at the end of the section compute a 
formula for the gluon distribution function valid to o? s . To compute the distribution 
function we need to sum over all possible color configurations. This color average 
involves the two dimensional gauge fields U. We start by recalling the properties of 
these gauge fields under our color averaging. 

3.1 Correlation functions involving the gauge transformations U 

According to equation (||) the gauge transformations U(xt) carry the information 
on the background field which enters the Green's function (|0|). These gauge trans- 
formations have the interesting property that the color average with the Gaussian 
weight (defined by equation @), of the combination 

UHxt)T a U{x t )U\y t )T a U{y t ) (18) 

can be written as Q 

<TrU\x t )r a U{x t )U\y t )T a U{y t )>= ^ ~ l) T {x t -y t ) (19) 

where N c is the number of colors and summation over repeated indices is implied. 

As pointed out in reference the average over the color sources yields the 
information about the ground state properties of the system. The average with the 
Gaussian weight is an artifact that simplifies the computation and we expect that 
as long as we resolve the nucleus on a transverse size much larger than the typical 
transverse quark separation, such an artifact is justified. 

The function T factorizes the dependence on the transverse coordinates xt, yt 
and is a function of their difference. Moreover, from equation (|l^) we see that T is 
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real and also 



T(0) = 1. (20) 
Denning the Fourier transform of T(xt) 

y(p t ) = J (fxte-^Tixt) (21) 
we have, together with (|20|), the sum rule 




The color charge at a given transverse location will be zero on average and the 
only way to generate a non-zero color charge will be by fluctuations. Equation (||) 
can be thought of as the generator of those fluctuations and thus the function 
F(xt,yt) represents the correlator of fluctuating fields at the transverse locations xt 
and y t . 

In momentum space, the function "f(pt) can be formally computed by expanding 
the exponential in (g) in powers of the coupling parameter a s \xjpt (weak coupling 
regime). This was done in reference || for scalars. For gluons, the result for 

a s p < Pt is 

7 (^) = (4vr) 2 ^iV c . (23) 
Pt 

We notice that the expansion is only necessary in order to analytically compute 
expressions of the form 

/ d 2 p t f( Pt M Pt ) , (24) 

with f(pt) a non-trivial function of pt- However, in principle, we can perform a 
numerical analysis to take into account the many possible different configurations 
of the external field contributing to expressions such as (^4|). This is equivalent 
to considering the effect to all orders in a s p/ P t of the different configurations of 
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the background field. For a quantitative discussion about the properties of the 
distribution function, we will restrict ourselves to the weak coupling regime for 
which j{pt) is given by equation (p3|). 



3.2 The distribution function 

With the above remarks in mind, we proceed to the computation of the gluon 
distribution function. We use the formula for the distribution function 



dN _ . 2k+ r dk- { dk' 

d 3 k ~ (2vr) 3 k+^k'+ 



/1£/^<W»' < 25 > 



where < Dff"(k, k') > is the small fluctuations propagator in momentum space, 
traced over the color and Lorentz indices and averaged over the external sources of 
color charge. This formula follows from computing < a}(p)a(p) > as an expectation 
value for the gluon Fock space distribution function in the ground state generated 
by the external valence charges H. 

In appendix I, we derive explicitly an expression for < Df^(k,k')>. Using this 
result, we obtain the following integral expression for the distribution function 
1 dN 2ik + 2 fdp + d 2 p t dk- 



irR 2 d?k (2tt) 3V c 'k+^k'+J (2vr) 4 



p\ — 2p + k~ — ie \ k~k + k' + 

1 1 



(^-(27r)V 2 )( Pt -fc, 
l(Pt ~ h) — 



p + — k + + ie p + — k' + — ie 
1 1 \ 



k + — iep + — k' + + ie 



u f 1 (n k t+Pt | (Pfh) 2 \ 

+ llPt l£ V k 2_ 2p +k--ie[ k-k+ + (k- 2 k+ 2 )J 

1 ' ' + _+ hrr^-T ^i— -77 1 1(26) 



\p + — k + — ieq + — k + — ie p + — k + + ie q + — k + + ie 

In the last term of this equation, we have taken the limit that k + — > k' + since this 
term has no singularity in that limit. 
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We now do the integral over k . We assume that k + > 0. When we do the 
integral, two classes of terms result. The first set of terms arise from the explicit 
k~ dependence in the above equation and are non-zero. The second set of terms 
arise from the q + in the last terms of the above equation. These terms result in an 
unrestricted integral over p + . One can show that all the singularities of the resulting 
integrand are on the same side of the p + integration contour in the complex p + plane. 
They therefore integrate to zero. (There is a possible ambiguity in the closing of 
contours associated with the contour at infinity, but this term does not have any 
contribution proportional to ln(l/x).) Therefore, we only get the contribution from 
the first term, which is only nonzero for p + < 0, 
1 dN , A 2 ^ 4fc+ f° dp + d 2 p t 



(N£ - 1 



TrR 2 d 3 k y c , (2ir) 3 J- 0O (2vr 



1 7( ^_^)-(2^)V 2 )(^-fe i ) 



p+(p+ - k+)(k+ + p+pi/kf) 
1 _ Pi I i + $\ + 2 ( p+ ^ 2 ( pt ' kt ^ 



(27) 



k+ \ k 2 t ) \k+ ) k? 
Now in this expression, we shall only be concerned with those terms which are 
proportional to ln(l/x). The terms not proportional to ln(l/x) are non-leading 
for small x. Moreover, we have found that within our approach, these terms are 
inherently ambiguous. This is due to the fact that the ln(l/x) terms can only arise 
by regulating the singularity in the above integral as p + — > oo. We do this by making 
the upper limit of integration, to be of the same order than the total momentum 
of a typical nucleon in the nucleus. Of course different regularization schemes will 
affect the non-leading terms in different ways. Presumably, the detailed longitudinal 
structure of the valence quark charge distribution must be known before these terms 
may be evaluated. 

After some straightforward algebra, we find 
( 1 dN \ 8(iV c 2 -l)l f d 2 p t r 



ttR 2 dxd 2 k t / (2vr) 4 x 

/ q 



J w? l{pt ~ h) [ l ~ -wr\ u 



(28) 
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where the subindex q in the left hand side of the above equation refers to the 
correction to the distribution function from the small quantum fluctuation field. 



Equation (28) is our main result. It is normalized so that the vacuum density is 
zero. This can be checked by setting U = 1 in the above equation. 
The terms above can be written in the form 



ttR 2 dxd 2 k t ) q \ 2vr 4 j x J (2tt) 2 ,yn ll \ p 2 k 2 



Now, shift the variable of integration pt — > pt+kt and expand ^(pt) in weak coupling. 
This restricts the lower limit of integration for the radial component of pt to be a s fi, 
which comes from the weak coupling expansion of 7. Thus the above expression 
becomes 

( I dN \ OA7/Ar2 ,,a 2 n 2 l 



8N C {N 2 - 1) 



ttR 2 dxd 2 k t J n CK c ; (vr) 2 x 

JO Ja 3 fi Pt [Pt + kf + 2p t k t cos 6\ \x 



where is the angle between pi and kt . The integral above can be performed exactly 
and the contribution to ([h]) from the ln(l/a;) terms becomes 



C(kt) = 2 ( ln( A) + I) . (32) 



with C(kt) given by 



Notice that the above expression means that the second term in the perturbative 
expansion of dN / dxd 2 kt in a s , develops the large factor ln(l/x) and that in the 
kinematical region of interest, the product a s ln(l/x) is not small. 

Furthermore, let us impose ordering in transverse momentum. This is the 
statement that the main contribution to the distribution function comes from the 
momentum region for which the emitted gluon's transverse momentum is larger 
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than that of the original one [10]. The effect is to restrict the integration interval 



for the radial component of pi in equation fl30|) which now runs between a s \x and 
kt- The reader can check that the above results in the modification of (j32^) which 
now reads like 

C {kt) =21nf-^-V (33) 

We can now include the contribution from the background Weizsacker-Williams 
field as given by equation (|i~3|). Thus finally, the perturbative expression for the 
gluon distribution function to second order in a s becomes 

Equation (|34| ) contains both ln(l/x) and ln(fct) corrections to the l/(xkf) dis- 
tribution and they represent the first order contributions to the perturbative ex- 
pansion for the distribution function. In the kinematical region of validity, these 
corrections are large. This signals that in order to properly account for the per- 
turbative corrections one has to devise a mechanism to isolate and sum up these 



leading contributions. Also notice that equation (28) is more general. In principle, 
it contains the information about the non-perturbative corrections as well. That 
information is in the function y(pt) an d it can be extracted by means of a Monte 
Carlo analysis for the whole kt domain. These issues will be treated in a future 
work. 

Diagrammatically, we can represent the background gluon field coupled to the 
external source (valence quarks), in momentum space, by means of figure la. The 
background field (wavy line) is by itself of order 1/g, according to equation (||) and 
the coupling to the external source (cross) can be considered to n-th order in the 
parameter g 2 fi/kt by means of the weak coupling expansion of equation (|8|). As an 
example, the Weizsacker-Williams distribution is obtained through the correlation 
of the background field taking the average over the source to first order (n=l) in 
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g 2 n/kt- This can be represented as in figure lb where the broken wavy line means 
that the momentum k is not integrated over. 

The gluon propagator in the presence of the background field can be computed 
perturbatively in the coupling constant g and the m-th order gluon propagator 
can be represented as in figure 2a. This is because the perturbative expansion of 
the gluon field involves its coupling to the background field through the covariant 
derivative and the background field acts as the gauge field. Notice that m has to be 
even since the gluon propagator is the correlator of two gluon fields and each time 
we couple the background field to the gluon field we introduce one power of g. In 
particular, the gluon propagator to second order in g can be represented as in figure 
2b. As suggested in this figure, the explicit dependence on a s of a quantity such as 
the gluon distribution function (which involves the gluon propagator) comes about 
only after performing our color average through the expansion in g 2 n/kt. This is 
because the coupling constant dependence of the background field and the order of 
the perturbative expansion offset each other. As shown in section 3, when computing 
the leading small x terms for the gluon distribution function, any term for which we 
can use the sum rule ( |22"|) will not exhibit an explicit coupling constant dependence. 
This becomes the criterion to decide that such terms are vacuum contributions. 

The gluon distribution function computed in section 3 can be represented by 
the diagram in figure 3, where we expanded the coupling with the external source 
to first order in g 2 ^i/kt- 

4 Loop corrections to the classical field 

Thus far, we have been concerned exclusively with the contribution of the small fluc- 
tuations propagator to the gluon distribution function. We have shown that this 
propagator induces large corrections proportional to a s ln(l/x) ln(fcf ) and a s ln(l/x) 
to the distribution function and have argued that the presence of these large log- 
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arithms signals the need to devise a method to sum them up to all orders in the 
perturbative regime. Before we do that we need to consider another contribution, 
to the same order, which comes from the corrections to the lowest order classical 
field induced by quantum fluctuations (see figure 4). This is apparant from equa- 
tion (|ll]) and equation ( |T^ ) where one sees that there is a contribution (A^ A^) p 
of the same order as ((A q A q )) p . 

In this section, we will compute the correction to the lowest order classical 
field induced by the quantum fluctuations. We will start by writing the total field 
A^ in terms of background (classical) and fluctuation (quantum) pieces allowing 
for the possibility that the background field may now be different from our lowest 
order classical (Weizsacker— Williams) solution. We will then write the equations of 
motion in terms of these new background and fluctuation fields keeping terms up 
to and including second order in the fluctuation fields. 

Our strategy will be to consider the expectation value of the equations of motion 
(in the path integral sense) and to relate the correlator of two quantum fields to 
the gluon propagator in equation (|T6|). We will show that only the + component of 
the equations of motion is modified and that the change could be thought of as the 
appearance of an induced current generated by the loop of fluctuation fields. We 
then proceed to explicitly compute this induced current and show that its effect is 
to renormalize the coupling constant g and the original background field. In other 
words, the modification induced by quantum fluctuations on the classical equations 
of motion can be cast into the standard expression for the renormalization of the 
coupling constant and the original background field to one loop in the light cone 
gauge. This result in itself is not surprising to QCD practitioners ( see for instance 
reference flTj). What is surprising is that this result persists to all orders in the 
effective coupling asH- 
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We start with the classical equations of motion 

W = 9J u a (35) 

and expand the full gluon field as 

A* = B» + b^ (36) 

where B^ is the background (classical) field, that is < > = B^ while is the 
fluctuation (quantum) field with < W > = 0. Keeping up to quadratic terms in 6 M , 
the + component of the equations of motion can be written as 

d-d-B- + (D t d-B\ = gj+ + g<J+> (37) 

where ja(%) = fabc < b\{x)d + b l c {x) >. Also, is the covariant derivative with B^ 
as the gauge potential. The corresponding expressions for the minus and transverse 
components of the equations of motion look like 

(D i D i B~) a - (Did-B\ + (D + d + B-) a + 
9fabc{ < b\d + b- > +(A < 6*6" >) bc + 

< bli&b-), >-< bi(D-b% > | = (38) 

(D 3 DiB l ) a - {Dj&B^a + (D + d + B% - (d-D i B~) a + 
(d-d-B\ + gf abc ^ < b\d + bi > +9_ < 67/6* > + 

(Dj < >) bc - < bl(D*V) c > + < y b {D' J b% > I = 0. (39) 

The expectation values of bilinear products of fields are related to the gluon prop- 
agator by the relation 

<bZ(x)bUy)>=-iG%(x,y). (40) 
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In the above, a, b,c. . . are color indices and /J,, v are Lorentz indices with i,j repre- 
senting the transverse components. 

The reader may verify that all the terms involving bilinear products of b^ in 
the minus and transverse components of the equations of motion either vanish by 
explicit computation, or, because they are symmetric in the color indices b and 
c, obviously do not contribute since they are always contracted with the totally 
antisymmetric structure constants / & c . In other words, the minus and transverse 
components of the classical equations of motion are not modified by the quantum 
fluctuations and the set of equations reduces to 

-d_d+B--(D i d+B i ) a = gjt+9<J+>, 
{D i D t B-) a -{D i d-B t ) a + (D+d + B-)a = 0, 
(D j DiB i ) a -(D j d i B*) a + (D + &+B i ) a -(d-D i B-) a = 0. (41) 

From now on we will concentrate only on the plus component of the equation of 
motion given by equation (^). It is clear that this equation is modified by the 
quantum fluctuations due to the presence of the induced current. In order to un- 
derstand this effect, we need to evaluate this term explicitly. For this purpose, we 
write it in the following way 

= f abc < bl(x)d + bi(x) >= if abc Km-^Gg(x,y). (42) 

Diagramatically, this term can be represented as in figure 4 where the wavy line is 
the background field and the spiral represents the loop of the fluctuation field. The 
loop is the vacuum polarization tensor and the component which contributes to the 
induced current and modifies the background field (which is purely transverse) is 
n + '. This allows us to represent the term d + G n = D + G u as H +l B l . 

We now proceed to compute the induced current explicitly. We will use our 
expression for the gluon propagator as given by equation ([!(]) . The first observation 
is that the terms in the Green's function with both x~ and y~ negative will be 
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symmetric in the color indices and will not contribute. Also, it can be shown that 
the terms with both x~ and y~ positive yield (after we implement the limit y — > x 
in a Lorentz covariant way) an infinite constant (independent of the transverse loop 
momentum) which vanishes upon dimensional regularization |18|| . 

However, the terms in the Green's function with opposite signs of x~ and y~ 
are a bit tricky. For these terms taking the partial derivative with respect to y~ 
followed by the limit y — > x is a very delicate operation and must be performed 
carefully. We find it more convenient to rewrite the terms with opposite signs of 
x~ and y~ in the Green's function in such a way as to avoid acting with d/dy~ on 
the terms 9(±y~). To do so, we will change the two dimensional integral over qt 
to a four dimensional integral over q. We can show that as a result, the product of 
theta functions of x~ and y~ will be replaced by theta functions of the light cone 
energy p~ . After some long but straightforward algebra we can rewrite the terms 
with opposite signs of x~ and y~ in the Green's function (which we call D |?) as 

d^, y) s (2.) / j^^n P - - 9 -)(p- + g- v4 )(g2 _, e) 

x{9(p-)e- l ^ x Tr[UHx t )T b U(x t )F c (p t - q t )\ x 4, + J^L-^^ - 1) 

[ |_ p p 

+ mj_r e iq+x- _\\ + QiPj(.4t-Pt) ,-ip+y- _ l\( e iq + x~ _ -p 

q~q + p~p + q~q + ' 



-e{-p-y^yTr[F\q t -p t )U\y t )T c U(y t )] x 



5 l3 + ^-(e^ + y- -1) 



+ PiPj_, e ip+x- _!)_)_ PiQjW ■ Pt) ,-iq+y- _ ]\( e *P + z- _ J> 

p-p+ p~p + q~q + 



(43) 



where F b {pt) is the Fourier transform of F b (xt) = U' (xt)r b U(xt). In the above 
expression, we are working with the color components of expression (|l6| ) (as opposed 
to the matrix notation), which is more suitable for the computation at hand. 
According to equation (^), we now have to compute 

3?(x) = ifabc^ x £zD b ^x,y). (44) 
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There are three distinct pieces to the computation corresponding to the different 
number of factors of p~ (q~) in the denominator of each term in expression fl43|) . 
For the rest of the section, we will outline the procedure for the computation of one 
of them, namely, the term proportional to 5{j and quote the result for the other 
two terms. In appendix II, we will show the detailed computation of some of the 
integrals necessary to fill in the intermediate steps. 

There are two important details to keep in mind while computing explicitly 
expression (f|4|): first, we have to implement the limiting procedure in a Lorentz 
covariant way. Second, when carrying out the integration over the relative and 
total (light cone) energies, we have to allow for a non-zero energy flow into the loop 
by not setting the total energy to zero, in spite of the structure of the integral which 
seems to require it to be so. The procedure is nothing but the well known point 
splitting method. 

In order to implement the limit y — > x in a Lorentz covariant way, we first 
transform the induced current to momentum space and then integrate over the 
relative (loop) momentum. We make use of the following identity. Let f(x,y) be a 
function of the four-dimensional variables x and y for which we want to compute 
the limit when x —* y. Fourier transform / to momentum space with respect to 
both x and y 

f(k,k') = J d A xd i y f(x,y)e ikx e ik ' y (45) 
and then make the change of variables 

s- k ~ k> 

where s and S can be thought of as the relative and total momenta respectively. 
Then 

f{s,S) = J d 4 xd 4 yf(.x,y)e ls(x - y) e iS(x+y) . (47) 
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Integrating over d 4 s/(27r) 4 will give e) 4 (x — y) which will set y = x upon integrating 
over y. Hence, 

' d*sf(s, S) = f d*x e i2Sx f(x, x) , (48) 



which is the expression for the Fourier transform of f{x, y) in the limit when y — > x 
as a function of 2S. With the above remarks in mind, we proceed to take the partial 
derivative with respect to y~ in equation ( |4"3| ) and then to Fourier transform with 
respect to x and y to the momentum variables k and k'. Let us look at the piece 
proportional to 5ij. Set i = j and then perform the p ± , and pt integrations to 
get 

9 r»«/-_ -A I n -\s,U~ , f d2 9t 



J d A xd%e ikx e ik 'y—Dl{x,y) {1) = {—2-K)8{k'~ + AT) | 



- k~)k' + 



(k'-)(k-) 

e(k'- 



2fe'~ A 2fc~ 

9(—k'~ 



(49) 



2fc~ 2fc' 

where the index 1 refers to our considering the first of the terms in equation 
namely, the term proportional to 5ij. Now, perform the change of variables ( f46| ) to 
write the above expression in terms of relative and total momentum and integrate 
over the relative momentum. The expression to evaluate becomes 

Tr[F b (s t + S t + q t )F c (S t - s t - q t )\ 



(2vr) 4 dy 
d A s d 2 q t {-2s-){S + -s 



(2vr) 4 (2vr) 2 (S~ - s~)(S- +s 

6(S 



( S + -s+- {S '-" t l 2 Z\ t )(S+ + s+- 



2(S"-s-) A" 1 a 2(S-+s~). 

0(-S- + 8-) 



(Q+ 4. s + _ {St+8t)*-ie - )( g+ _ + _ (£tz!f) 



(50) 



2(S-+S") A w 1 2(5- -s-)- 
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We will perform the transverse momentum integrations last. For the moment, 
let us concentrate on the s + and s~ integrals. It is more convenient to do the 
s + integral first since this can be done by contour integration. The s + dependent 
integral of the above expression, denoted by /, is 



This integral has a logarithmically divergent piece which can be isolated by adding 
and subtracting the term (St — St) 2 /2(S~ — s~) to the numerator. This divergent 
piece can be shown to give a constant independent of the transverse loop momentum 
and therefore it vanishes upon dimensional regularization in the transverse direction. 
The remaining piece reads as 



Let us investigate the above expression in some detail. The integrand has two poles 
in the complex s + plane. Their location depends on the signs of (S~ — s~) and 
(S~ + s~). If the two poles are on the same side of the real axis, then we can close 
the contour of integration on the other side of the real axis and the integral vanishes. 
So in order to get a non-vanishing result, the two poles must be on opposite sides 
of the real axis. Recall that 2S is the total or external momentum flowing into 
the loop and thus it has to be kept fixed (and finite for a non-trivially zero loop 
integral) while working in momentum space. Thus for a given sign of S~ only one 
of the two terms in the integral (|5^) above contributes. This can be seen as follows: 





(51) 




(52) 
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the first term in (52) is non-zero only if the two conditions 

S~ - s~ > 0, S~ + s~ > (53) 
are satisfied, whereas the second term is non-vanishing only if 

S~ - s- < 0, S~ + 8~ < 0. (54) 



First, take S~ positive. Then only the first condition (|53| ) gives overlapping in- 
tervals for s~ namely S~ > s~ and s~ > — S~ or — S~ < s~ < S~ whereas the 
second condition (54) does not. Therefore the second term in equation (f)2) can be 



disregarded and only the first one is non- vanishing. The opposite is true for S~ 
negative in which case only the second of the terms in equation ( |52|) contributes. 

The integral ( [50D however turns out to be independent of the sign of S~ (as we 
shall describe below). The reason is the scaling property of the integral over s~ . This 
can be understood by recalling that after all, the overall expression, equation (|50|), 
is explicitly proportional to 5{2S~) and any term proportional to S~ can be thrown 
away after scaling the integration variable s~ by S~ . With these remarks in mind, 
let us continue working with a definite sign of S~ , say, S~ > 0. Performing the 
integral (|5^) we get 

-iO(S- - s-)9(S- + s-)(S~ + s-)(S t - s t ) 2 

s-[as+s-(i- + -% 2 (i-0- (s t - St ) 2 (i + 0] 1 j 

in terms of which, equation ( |50| ) is 



(2vr)4 dy- 
(-2i)8(2S-) J ^^LTr[F a (s t + S t + qt )F b (S t -s t - qt )](S t -s t ) 2 

s~ds~ 



(56) 



-5- (S--s-)[4S+(S- + s-)(S--s-)-q 2 (S--s-)-(S t -s t ) 2 (S-+s-)} ' 

Let us now look at the integration over s~ . We scale s~ by S~ , that is we define the 
variable £ = s~ /S~. We notice that after scaling we can make use of the explicit 
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factor 5(2S~) in equation (|50j) and we can safely throw away any term which is 
still explicitly proportional to S~. Thus, the term [4S + S'~(1 — £)(1 + £)] in the 
denominator of the above drops. As a result, the overall expression ( |50j becomes 



/ 



^ S j4 ^4 is(x-y) iS(x+y) ® n ii i \ 

- { d xd ye v y 'e y y '-^D bc (x,y)^) 



(2tt)4 dy- 
(2i)5(2S-) J ^L^SLTr[F a (s t + S t + q t )F b (S t - s t - q t )](S t - s t ) 



(57) 



-i (i-OW(i-0 + (5t-*t) a (i + 01' 

To proceed further, it is convenient to shift qt — > qt~ st- Then the arguments of 
F' s in the trace of the above expression become independent of st and can be taken 
outside the st integral which will be evaluated next. The st integral is a formally 
divergent integral and must be regulated. This is done in appendix II using the 
dimensional regularization method. We find that 

fd^st (St - s t ) 2 _ r(-a;) 2 

which brings expression ([57]) to read like 



(2tt) 4 y Qy-"**-,*^!) v-v 16?r 

x / ^Tr[F a {S t + q t )F b {S t -q t )]{S t -q t ) 2 fedi. (59) 

The integration over £ can now be done easily. It just gives a factor of 2/3 and 
the remaining transverse momentum integral can be computed by using the explicit 
form of F in terms of the gauge transforms U. We also show this in appendix II 
where we find that 

(hqt ;Tr[F b (S t + q t )F c (S t - q t )](S t - q t f = £ f bcd p d (2S t ) (60) 



(27r) c v w t ii/jv^i — 2 



with 



p d (2S t ) = | tfxttP^piixt) (61) 
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being the Fourier transform of the charge density with respect to 2St- Putting 
everything together, we get finally the result that the expression for the Fourier 
transform of the induced current coming from the term proportional to Sij as a 
function of 2S is 

3+(2S) (1) = if abc J ^d*xd%e^-yh^ + y^Dt c (x,y) {1) 



8tt 



g 2 p d (2S t )5(2S-). (62) 



Above, we have written only the divergent part of the gamma function when u — > 1 
and have used that for SU(3) the product fabcfbcd = CA^ad with = 3. The 
remaining three terms in equation (E3) can be evaluated in a similar fashion. Here 
we just quote the result 

ia + (2S)(2) = T -^g 2 ~ Pd {2S t )5{2S-) 
ia + (2S) (3) = T -^g 2 ~ Pd {2S t )5{2S-) 

j+(2S) (4) = 4^-^ 9 2 ^(2^)<5(25-). (63) 

The final result for the Fourier transform of the induced current is obtained by 
adding the four terms given by equations (62) and d63| ) and it becomes 

j+(p) = g 2 T(-u;)^p a (p t )5(p-) , (64) 

where we have renamed 2S — > p. Thus the term that modifies the plus component 
of the equations of motion in momentum space just becomes gja(p)- 

We proceed to argue that by rewriting this result in terms of an expression 
involving the components of the polarization operator, we can absorb the effect of 
the loop corrections on the equations of motion, into the renormalization of the 
coupling constant. Let us first recall that according to equations (|j) and (||), which 
are the classical solutions to the equations of motion, the expression for the zeroth 
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order background field in momentum space, can be written as 
Af )(p) = f d^xe^^ix) 



= ^)9{j^)p b {pt)6( P -). (65) 
Therefore, notice that gj^ip) can be written as 

9jt(p) = Kt(p)4 0) (p) (66) 

with n^(p) given by 



fsr(-u) 



™ =9pV ll^J^ (67) 

which is the standard expression for the +i components of the polarization operator 



in light cone gauge l£j. The fact that we recover this well known result is truly 
remarkable and is one indicator of the success of our formalism. 

We can now take an ansatz for the formal solution of the system of equations ( |4l| ) 
to be 

B~{x) = 0, 
B\x) = 9(x-)<j R (x t ), 
t-o? r {x) = —UlxtytWixt). (68) 

gn 

where gn is the renormalized coupling constant whose expression is obtained 
through the computation of U^ u and will be given explicitly by 

g\ = Z z g\ Z S = ( 1 + ~ ) ■ (69) 



This in turn means, according to (|68D, that the field B gets renormalized by the 
inverse of the constant that renormalizes g: 

B i {x) = Z- 1 l 2 A i ^{x). (70) 
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The above exercise has taught us the important lesson that the modifications to 
the background field introduced by the quantum fluctuations do not induce extra 



terms in the expression for the distribution function (28). Further, their effect 
can be included by replacing the coupling constant g by the renormalized coupling 
constant gR. 



5 Summary 

We have presented in this paper an expression for the quantum corrections to the 
Weizsacker-Williams gluon distribution at small x valid to all orders in the param- 
eter a s fi. We used this expression to compute explicitly the leading ln(l/x) and 
ln(/ci) terms in the momentum regime a s fi <C kt- We have shown that the pertur- 
bative approach introduces a series expansion parameter a s \n{l/x) which is large 
and thus, forces us to devise a method to sum up the leading contributions to all 
orders in that expansion parameter. Nevertheless, the present result already signals 
that at small x values the gluon distribution function will be modified significantly 
from the l/(xk?) behavior. 

We wish to emphasize that our central result, equation (|28"|), contains in prin- 
ciple the information about the quantum correction to the classical distribution to 
all orders in the parameter a s /i. 

We have found that the only effect the quantum corrections have on the classical 
background field can be absorbed into the renormalization of the field and the 
running of the coupling constant. 

We have not addressed the issue of summing up the perturbative series in this 
paper. In the weak coupling limit, this is equivalent to solving an integral equation 
for virtual corrections to the gluon propagator. Another issue we would like to 
address is whether we can relax some of the constraints in the model. In particular, 
we need to pay attention to the restriction set by the necessity of having a large 
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(perhaps too large) A nucleus in order to compare our predictions to experimental 
data. These issues will be addressed in a future work. 

Appendix I: The propagator in momentum space 

In this appendix, we shall Fourier transform the small fluctuations propagator to 
momentum space. This result will be useful for computing the distribution function 
and may also be useful for other computations. We will first consider the propagator 
without averaging over all possible color orientations of the source fields. 

We shall derive two representations of the propagator. One representation will 
include an integration over the variable p + . This representation will be useful for 
computing distribution functions. We will also present a second representation 
where we have performed this integration over p + . Finally, we will present a result 
for the propagator after averaging over the colors of the external field. 

We recall that in coordinate space the propagator has the form 




) 





1) 



QiPjpt ■ qt 



( e -iv + y _ i)( e *<7 + * _ i) 



(71) 



(p p+)(p q+) 
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2 2 

where q + = p + + \ Pt and 



2 P - 

a/3 



F^(x t ) = (u(x t )r a U^x t )Y . (72) 

We now wish to define the Fourier transformed Green's function 

Gf' a 'P' (k,k') = j d A xd A ye- lkx+lk 'yGf' a 'P' (x,y) , (73) 

which we can divide into four pieces as 

Gf' a 'P'(k,k') = G a ^\k,k') + G a + ^'(k,k') 

+ Gl% P \k,k') + G a _% P \k,k') (74) 

In this equation, the first index ± refers to the index of 9(±x~) and the second to 
0(±y~) in the definition of the coordinate space Green's function. 
Let us first consider the component 

G°?j?/(k,k>) = - J ^iAe*- t ' I - i ( p - t '^(-i^(-f)r5;'^ 

,ip + (x~-y-) _ e -ip + y~ _ e «p+:r^ j _ 



x <^ Sa + 



in 

P V 



We can perform the integrations over p , pt, x and y to obtain, 

G a _%? {k,k') = -(2*)5{k- - k'~)(27T) 2 5^(k t - k' t )rfr^'A^(k,k') (76) 
where 

A (k, k - Jdx dy — 

x f^Sij + ■^ T (2e ip+{x '- y ' ) - e- ip+y ~ - e ip+x ~)^ (77) 

Now let us do the x~ and y~ integrations. In performing these integrations, ie 
factors will appear which will guarantee the convergence of the integrals at infinity. 
We find 



A (kk >) = f d P + 1 I 

v ' ' J 2ir k 2 t - 2p+k- -ie p+ -k+ + ie p+ - k'+ - ie 

* { 5 ^k^k^ i2P+ - k+ - k ' +) } 



(78) 
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Finally, we can perform the integration over p + to find 



k'+ -k+ + ie k' 2 - ie 

+ 



-{k' + - fc+)) + k,k' -> -k',-k^j (79) 



k~(k+ -ie){k'+ + ie) 

Now, to fully define this Green's function, we must specify the nature of the 
singularity at k~ = 0. In the last equation, we would like to use the Leibbrandt- 
Mandelstam prescription on l/(k~ + ie/k + ) whenever we have the combination 
l/k~k + and l/(k~ + ie/k' + ) whenever we have the combination l/k~k' + . We 
would like to go further however and define l/k + k~ as l/(k~k + + ie). This can be 
done as follows: we use l/{k + — ie) = 2iri5(k + ) + l/(k + + ie) whenever we have the 
constraint that k~ > 0, and a similar modification when k~ < for k' + . 

On the other hand, it is more difficult to implement the Leibbrandt-Mandelstam 
prescription in the expression which involves the integral over p + . However, we will 
only use this result when both k + and k' + have the same sign and are non-zero. In 
this case the Leibbrandt-Mandelstam prescription is unambiguous. 

Our results for the two representations are therefore 

G^f(M') = -2 m 5(k--k-')(2^)H^(k t -k^Tfr^' I ^— 



k'+ - k+ + ie 



»(*-) (« tf + F ^ S :(^-* + ) 



+ k,k' -► -k',-k\ +2m5{k + )6{k-^^ L 



k~k'+ 

+ k,k' -> -k',-k\ (80) 

Here the k + , k' + and k~ singularities are treated using the Leibbrandt-Mandelstam 
prescription. 

For k + and k' + both non-zero and of the same sign, we have the representation 
G a _% P \k,k') = -2^{k--k-'){2^^\k t -k' t )rfr^' 
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/ 



dp + 1 1 1 

2ir kf — 2p + k~ — ie p + — k + + ie p + — k' + — ie 

where the k~ singularity is treated using the Leibbrandt-Mandelstam prescription. 

The evaluation of the remaining contributions to the Green's functions can be 
done by the same methods as above. There is nothing really new in the analysis 
except that it is longer and more involved. The subtlety is in the treatment of the 
singularities in the k^ and k' + variables. This has been discussed above and is 
treated as such. The results are 

G a + % P \k,k') = -2m6(k--k-') J ^L F ^( Pt -k t )F^'(k' t -p t ) 



1 



k+ - k'+ + ie 



(2vr) 2 



~\~ /tj k ^ k k 



+ 



k-k+k'+ y ') p\ - 2k~k+ - ie 

27ri5(k' + )0(k-)^- ! 



' k k+ pi — 2k~k + — ie 
+ (82) 

For k + and k' + both of the same sign and non-zero, the above is equivalent to 

cSf (*> k ') = ~^ k ~ - k ") J ^N*^ ~ ~ ft) 

1 1 1 



p\ — 2p + k — ie p + — k + — ie p + — k' + + ie 

* (^ + *w< 2 " + -* + -*' + >)- (83) 

We finally also obtain an expression for It turns out that in this expres- 

sion, no restrictions on the values of k + and k' + are needed to get the singularities 
in 1/k^ or l/k' + into the Leibbrandt-Mandelstam form. The results are 

G%%P\k,k') = 2m5(k- - k -')2k-9{k-)^— J ^F^(p t - k t ) 

X K {kt Pt) pi-2k-k+-ie\ dl3 k~k'+ k-k+ 
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Pik'jpt • k' 1 

+ {k-k+)(k-k'+)j ' (84) 
We also have the equivalent integral representation where 

Gflif (k, k>) = 2,5{k~ - k-') J dj ^F^( Pt - k t )F^'(k' t - Pt ) 

1 1 1 



kf — 2k p + — ie q + — k + — ie p + — k' + — ie 

\° 13 k-k'+ k-k+ + (k-k+)(k-k'+)j- { ' 



In this equation, q + = p + -\ — 
The expression for G_ + is 

&ff/(k, k!) = G a _%P\-k\ -k) (86) 

We now want to convert these expressions for the propagator from the matrix 
basis to the component basis. To do this, we make the transformation 

{U{x)T c U ] {x)){U{y)T c U ] {y)) -► A{Tr T a U{x)r c U ] {x)){Tr T b U{y)T c U ] {y)) (87) 

We then use the identity 

T a(5 T a>f3> = ^ (<W"^'/3 ~ Jj-S a/3 S a > p^j , (88) 

which results in following transformation for the definition of the propagator 

F^(x)Fl? p/ (y)^2TrF a (x)F b (y) (89) 

We now want to proceed to derive formulae for the propagator which has been 
averaged over all values of the color charges of the valence quarks. To do this we 
define 

< Tr U ] {x)T a U{x)U\y)T h U{y) >= ^T ab (x - y) . (90) 
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Notice that 

T ab (0) = S ab 

We now define 

< G$(k,k') >= (2ir)<5(fc _ - k f -)(2Tr) 2 5(k t - k' t )D? b {k,k') . 
As before, we can write Dff as 

D% = D± l3 {k, k!) + Df +ij (k, k') + Df^k, k') + Dl b +ij (k, k') . 
Using the above substitutions we find that 

-,ab rab [ dp + I 



(91) 



(92) 



(93) 



D° 



-u 



X 



5 a + 



2tt kf — 2p + k~ — it p + — k + + it p + — k' + — it 



A/a k 



i w 3 



-(2p+ - k + - k'+) 



(94) 



'■' ' k /.• //• 

which is valid for k + , k' + nonzero and sign{k + ) = sign{k l+ ) . The expression with 
the p + integral completed is 



D 



ab 



+ 2md(k + ^^- 



1 



-v 



For -C++ we obtain 

dp + d 2 pt 



k'+ -k+ + it 
~\~ kj k ^ k j k / . 



5ij + L-k+k>+ 



kikj (k' + - k + ) 



(95) 



n ab . . 



(2vr) 3 
1 



r ab (p t - h 
i 



p\ — 2p + k — it 



PiPj 



(2p + - k+ - k' + ) 



,(96) 



p + — k + — it p + — k ,+ + it 
which is valid for k + , k l+ nonzero and sign(k + ) = sign(k f+ ) . After completing the 
p + integral, 

1 



p 2 — 2k~k + — it 



1 



k+ - k'+ 



PiPj 



k~k+k'+ 



(k+-k ,+ )] + 27Ti5(k' + )^L 
) k k + 



~\~ k^ k ^ k , k 



(97) 
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We also have 



D ab ■■ 



dp + d 2 p t b 



r ab ( Pt - h] 



(2vr) 3 ' '" LJ k 2 - 2k~p^ 



ie q 



+ 



PikjPt ■ h 



k~k'+ k-k+ (k-k+)(k-k'+) 



k + — ie p + — k' + + ie 
(98) 



where, again, q + = p + + 



+ , Pt- k t 



2k- 



. Doing the integral over p + gives 



D ab ■■ 



(2*) s 
PiPj 



1 



t} p\ -2k~k+ -ie 
PikjPt • h 



13 k~k'+ k-k+ {k-k+){k-k'+) 



Finally, the expression for D |_ is given by 



(99) 



Dl b +ij (k,k') = D b + a _ ji (-k',-k). 



(100) 



Appendix II 

In this appendix, we will explicitly evaluate two of the integrals we use in the 
computation of the one loop corrections to the classical background field. The first 
integral we will evaluate is 

1 J (2vr)2 [(q t - St ) 2 (l - + (St ~ s t ) 2 (l + 0] ' [ ' 

Let us first shift st — > st + St and write the above as 

d 2 s t 

(2vr)2 l( St + S t - qt ni-0 + ,f(\ +0] 
The denominator in the integral can be written as 



2s 2 + (2s t ■ (S t - q t ) + (S t - qt?)(l ~ 0- (103) 



Performing the change of variables 

, (l-0 ,o 
v t = s t -\ 2 — ( S t-Qt) 

d 2 v t = d 2 s t 



v 2 = s 2 t+Sf(St- q t)(l-0+ { -^^(St-qt) 2 (104) 
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Eq. ( |103| ) can be re-written as 



^[v 2 t+^—^{S t - qt f] (105) 



and the integral ( |102| ) 



''J ^) 2 2[v 2 + ^(S t - qt ) 2 ] (Ubj 
Written in this form, we can drop the linear term in vt and the contributing terms 
to flloel ) will be 

d 2 v t v? 



(2vr) 2 2 [v 2 + - Qt ) 2 ] 

r £ Vt (i - pyA(s t - qt f 

The first and second terms in expression ( |107| ) are quadratically and logarithmically 
divergent. To regulate the divergences we compute them by dimensional regular- 
ization. We thus write ( |107|) as 



_ 1 ( 2tt 2 \ J r°° vj 

2(2ny{T(d/2))\Jo Vt [vf + ^fl{S t - qt f] 



4+i 



+ Jo ** y t + a^i {St - qt?] I {m 

To evaluate the integrals in the above expression we use the well known formula 

u?du r(i(l + /3))r(a-i(l + /3)) 



Jo (u 2 + C 2 Y " 2(C 2 ) a -( 1 +l 3 )/ 2 F(a) 
by means of which (|108|) becomes 



h = (^)(^) 2 2vr d / 2 r(- ( i/2)( T ^-)[(l - f)(S t - qtf/Af 2 . (110) 

We are interested in the divergent part of this expression when d — * 2 thus for that 
we take d = 2 everywhere except in the argument of T. We obtain finally 

J 1 = 1^2^(1 -ft) 2 , (in) 
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Next, we want to compute the following integral 

h=J ■^Tr[F a (S t + qt )F b (S t - qt )][(S t - qt ) 2 -(S t + q t ) 2 ]. (112) 

For this purpose let us write F in its explicit form in terms of U. 

F a ( Pt ) = J d 2 x t e ip * Xt U\x t )r a U(xt), (113) 

and thus the integrand in ( |112| ) can be written as 

' d 2 x t d 2 y t e i ^ + ^ Xt e i ^-^ yt [(S t - qt f - (S t + qt f] 

Tr[U\x t )T a U(x t )U\y t )T b U{y t )}. (114) 

The factors (St±<jt) 2 are to be interpreted as derivatives acting on the corresponding 
exponential. Integrating by parts, ignoring the surface terms and with the help of 



Eq. (114) we can write Eq. (|112j) as 



I 2 = - J ^ d ^ xtd ^(s t +^ te i(s t - qt )y t 

x {Tr[U\x t )T a U(x t )V 2 (U\y t )r b U(y t ))] 
- Tr[V 2 (W(x t )T a U(x t ))uHy t )r b U(y t )]}. (115) 
Performing the qt and yt integrals, Eq. ( |112| ) becomes 

h = - [ d 2 xte 2iS ^{Tr[U\xt)T a U{x t )V 2 {U\x t )T b U{xt))] 



- Tr[V 2 (uHx t )T a U(x t ))U^x t )T b U(x t )]y (116) 

Since the gauge transformations U and the charge density are related by Eq. (g), 
one can prove the following identity 

Tr[[/WV 2 (£/W) - V 2 (C/W)£/W] = _^/o6c pc(a . t) _ fbac pc{xt) y n7) 

Using the antisymmety of f abc and plugging the above back into Eq. ( |116j ) we finally 
obtain the result 

h = g 2 f abc J d 2 x t e 2tStXt Pc (x t ) , (118) 
which is the Fourier transform of the charge density with respect to 2 St. 
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Figure captions 

Figure 1: a) Coupling of the classical background field to the external source. The 
wavy line represents the background field which is of order O (^). The external 
source is shown by a cross and to the lowest order in the weak coupling regime it 

2 

is O (^-)- b) Correlation of two classical background fields where the broken wavy 
line means that the momentum kt is not integrated over. 

Figure 2: a) The perturbative expansion of the gluon propagator in the presence 
of the classical background field in terms of the coupling constant g. b) The gluon 
propagator expanded to the second order in g. 

Figure 3: The correction to the Weizsacker- Williams distribution function to the 
lowest order in the weak coupling regime. 

Figure 4: Modification of the classical background field due to quantum fluctuations. 
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